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Abstract 

We compute the gamma-ray output of axion-mediated dark matter and derive 
the corresponding constraints set by recent data. In such scenarios the dark matter 
candidate is a Dirac fermion that pair-annihilates into axions and/or scalars. Provided 
that the axion decays (at least partly) into photons, these models naturally give rise to 
a box-shaped gamma-ray spectrum that may present two distinct phenomenological 
behaviours: a narrow box, resembling a line at half the dark matter mass, or a wide 
box, spanning an extensive energy range up to the dark matter mass. Remarkably, we 
find that in both cases a sizable gamma-ray flux is predicted for a thermal relic without 
fine-tuning the model parameters nor invoking boost factors. This large output is in 
line with recent Fermi-LAT observations towards the galactic centre region and is on 
the verge of being excluded. We then make use of the Fermi-LAT and H.E.S.S. data to 
derive robust, model-independent upper limits on the dark matter annihilation cross 
section for the narrow and wide box scenarios. H.E.S.S. constraints, in particular, turn 
out to match the ones from Fermi-LAT at hundreds of GeV and extend to multi-TeV 
masses. Future Cerenkov telescopes will likely probe gamma-ray boxes from thermal 
dark matter relics in the whole multi-TeV range, a region hardly accessible to direct 
detection, collider searches and other indirect detection strategies. 



1 Introduction 



There is mounting evidence for the existence of dark matter (DM) in our Universe through its 
gravitational interactions with ordinary matter (see for instance [III2] for reviews). However, 
there is currently no unequivocal indication for non-gravitational interactions neither in 
direct, indirect nor collider searches. Good prospects exists, though, if the dark matter 
particle is a thermal relic with a mass of the order of the electroweak symmetry breaking 
scale In this case, the interaction rate with nuclei, the annihilation rate in galaxies 
or the production rate at colhders could be large enough to allow detection, provided the 
theoretical and experimental backgrounds are sufficiently well-understood and suppressed. 

A promising strategy to detect dark matter particles consists in the search for gamma- 
ray spectral features produced in dark matter annihilations. Such features could stand out 
over the featureless energy spectrum of the diffuse gamma-ray background or the known 
astrophysical sources even for moderately small annihilation rates. Up to now three sharp 
spectral features have been identified: gamma-ray lines [H El E] , internal bremsstrahlung [TJ 
[HI [9] and gamma-ray boxes [10] . Gamma-ray lines appear in the annihilation or the decay of 
dark matter particles into two daughter particles, one of which is a photon, hence producing 
a spectrum with a monoenergetic line. Internal bremsstrahlung occurs, for example, in 
annihilations of Majorana or scalar dark matter particles into a photon plus a relativistic 
fermion-antifermion pair, giving rise to an energy spectrum with a prominent feature close 
to the endpoint. Lastly, a gamma-ray box arises in the annihilation or decay of dark matter 
particles into a pair of intermediate scalars, which in turn decay into two particles, one 
of them a photon. Then, the photon, which is monoenergetic in the rest frame of the 
intermediate particle, is boosted in the galactic frame and acquires an energy that depends 
on the angle between the momentum of the photon and the momentum of the intermediate 
scalar, resulting in a box-shaped spectrum. 

Several searches for these gamma-ray features have been conducted in recent years. 
Gamma-ray lines have been searched for in the galactic centre [TTl [T21 US] , in the isotropic 
diffuse gamma-ray background [SI US] and in dwarf galaxies [12], probing at present the 
range of dark matter masses between 1 GeV and 20 TeV. Besides, signatures from in- 
ternal bremsstrahlung have been searched for in the galactic centre [Ml [13] and in the 
isotropic diffuse gamma-ray background [13], covering the mass range 40 GeV — 300 GeV 
and 500 GeV — 20 TeV. Finally, gamma-ray boxes have been searched for in the galactic 
centre and halo regions [lU] for dark matter masses in the range 5 GeV — 800 GeV. No sign 
of dark matter annihilations has been unequivocally found, although an intriguing hint for 
a sharp feature at ~ 130 GeV has been reported in [161 HZl [IE]- It is at the moment unclear 
whether this feature is indeed the result of dark matter annihilations or, on the contrary, 
is an instrumental artifact or a statistical fluke. Future observatories have, however, good 
prospects to ellucidate whether this gamma-ray excess is genuine or not |19j . 

From the theoretical side, previous works have discussed the phenomenology of particle 
physics models which generate gamma-ray lines [20l[2ni22l[251[2il[25l[26l[27l[28ll29], 
signals from internal bremsstrahlung [301 1211 [321 [SI [M] or gamma-ray boxes [TOl [291 ES] . 
In this paper we shall focus on the axion-mediated dark matter model recently proposed 
in Refs. [2S1 [2Z1 [3B] and discuss a speciflc realisation that generates a box-shaped gamma- 
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ray spectrum. This spectrum falls naturally into one of two possible categories: a narrow 
box or a wide box. In both cases thermal relics produce gamma-ray outputs at the reach 
of current experiments. The search for gamma-ray boxes is thus well-motivated from the 
theoretical viewpoint and, on that respect, the next generation of very high-energy gamma- 
ray telescopes (such as CTA [32]) is eagerly awaited. 

The paper is organised as follows. In Section |2] we introduce the main features of the 
axion-mediated dark matter model and in Section [3] all the relevant tree-level annihilation 
cross sections are explicitly presented. Then, in Section |4] we compute the induced box- 
shaped gamma-ray spectra and derive model-independent constraints on the narrow and 
wide box scenarios using Fermi-LAT [ID] and H.E.S.S. [H] data. Our final remarks are 
drawn in Section IH 



2 Axion-mediated dark matter 

We consider the model presented in Refs. [361 EH EH], where the Standard Model (SM) 
is extended by a Dirac dark matter particle x ^ complex scalar S charged under a 
global [/(l)pQ symmetry. Similar models were studied in Refs. [121 HSj in the context 
of a dark matter explanation for the PAMELA positron excess. The Peccei-Quinn (PQ) 
transformations of the dark matter particle and the mediator are, respectively, x ~^ e*'^^"x 
and S — )■ e~^*°S'. We further assume that S acquires a vacuum expectation value such that 
the complex scalar field can be cast as S* = (f ^ + s + ia) / v^, where Vs = V2 (S) is the axion 
decay constant. The effective Lagrangian of the model then reads 

^ = ^xr^^X + \d^S\^ - V{H, S) + Ant + Csm , (1) 



where £sm is the Standard Model Lagrangian, which includes the Higgs potential Vsm 
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Xh\H\^ + m%\H\'^, and 



V{H,S) = \s\S\' + 2\Hs\Sm\' + ml\S\'-{^m'^S' + c.c.), (2) 
Ant = -XASxPLX + S*xPRX)+Y.I^''^i^'^^'"' (3) 



1=1,2 ^ 



with F^i, = \e^ypaFP'^ . The constant parameters q depend on the charges of the heavy 
fermions in the anomaly loop diagrams and are explicitly given in Refs. [371 EH]- It is 
possible to show that the Lagrangian contains a mass matrix for the CP-even scalars which 
can be diagonalised by the field rotation 



s = cos 6 s + sin 6 h, h = — sin 6 s + cos 6 h , (4) 
where the mixing angle is given by 



, „2 2XhsVsV . . 



Vs and V being the vacuum expectation values of the singlet and the Higgs doublet. 
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The axion has a mass nia = m'g which stems from the mass term ^m'gS'^ + c.c. This 
mass term can originate from the higher dimensional operator $^5*^ + c.c, where $ is 
a high-scale Peccei-Quinn breaking field, with ($) = fa ^ Vg- The high-scale PQ breaking 
scale is constrained to be 10^ GeV < fa < 10^^ GeV from bounds of supernova cooling and 
relic density on the invisible axion coming from $, respectively. Thus, the weak-scale axion 
mass rria ranges between 1 GeV and 10^ GeV, for g ~ 0{1). 

The interaction Lagrangian in Eq. (|3| gives rise to a dark matter mass and to dark 
matter interactions with the CP-even scalar s and the CP-odd scalar a which, assuming for 
simplicity that the coupling A^^ is real, are given by 

^DM = -m^XX- ^K^XX- ^K^^X^^X, (6) 

where rriy^ = X^Vg/ a/2; the interaction term with the mass eigenstates 5, h can be easily 
obtained from Eq. Q. 

Lastly, the effective axion interactions can be rewritten in terms of physical electroweak 
gauge bosons: 

'^aVV = CV^Vjaf^iiupuFy^ F^" , (7) 
where the only non-vanishing terms are 



(ciOi cos 9w + C2012 sin Q-w) 



s 



cz-y = 77 (C2a2 - cioi) sin(26'H/^ 



czz = 77 (02^2 cos^ 6w + ciai sin^ 6w) , 



_ C2a2 

CWW — a • 1,0; 

With these interactions, the dark matter particle is absolutely stable due to the residual 
Z2 symmetry which emerges after the Peccei-Quinn symmetry is spontaneously broken by 
the singlet vacuum expectation value. However, the CP-even singlet scalars s, h can decay 
into SM particles due to the mixing with the Higgs boson. Additionally, the heavier eigen- 
state s may decay into a pair of the lighter eigenstates h, and both eigenstates s, h may 
decay into a pair of dark matter particles if kinematically allowed. The decay widths for all 
these processes can be found in Ref . [30] ■ We also note that the scalar eigenstates may also 
decay into a pair of axions if kinematically allowed. 

The axion, on the other hand, can decay into two electroweak gauge bosons with rates 
given by 
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Figure 1: Branching ratios of the axion decay as a function of the axion mass nia for a ratio of 
anomaly coefficients ci/c2 = 1/3, 1, 3, when the axion couples to vector-like doublet {Y = ±1/2) + 
triplet {Y = 0), vector-like doublet {Y = ±1/2), and vector-like doublet {¥ = ±1/2) ± vector-like 
singlet {Y = ±1) fermions, respectively. Solid black, dashed grey, dashed blue and dot-dashed red 
lines correspond to the WW, ZZ, 77, Z7 decay modes, respectively. 
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Decays into Standard Model fermions are not allowed at tree level. 

The branching ratios for the axion decay as a function of the axion mass are shown in 
Fig. [T]for various choices of Ci/c2. It is apparent from the figure and from Eqs. (8][9) that, 
for generic choices of parameters, all decay channels have sizable branching ratios provided 
they are kinematically open. More specifically, the decay channels a — )■ 77 and a — )• Z7 can 
have large branching ratios and even be the dominant decay channels. This property of the 
model will have important implications for indirect searches, as will be discussed in detail 
in Section HI 



3 Tree-level annihilation cross sections 

At tree level there are three possible dark matter t-channel annihilations, XX a s, a a, s s, 
which easily dominate over the one loop induced s-channel annihilations into gauge bosons 
XX V1V2 unless ~ 2m-^, namely when the rate for the latter processes is resonantly 
enhanced. In this section we will derive the velocity-averaged annihilation cross sections for 
these processes assuming for simplicity that the mixing between the two CP-even scalars 
is zero; the full expressions including the Higgs-singlet mixing can be easily derived from 
the formulas below using the mixing angle defined in Eq. ^ and the mass eigenvalues 
corresponding to s and h. 

The dark matter annihilation cross sections times relative velocity (without thermal 
average) for all the t-channel processes are given by 

1 / (ma±m,)2\i/2/. (m^ -ms)2\i/2 1 r+i 
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(11) 
(12) 



where Eqm = 2m-)^{l — f^/4)^^/^ and 9 is the angle between the direction of the incoming 
dark matter particle and the outgoing axion or scalar, both expressed in the center of mass 
frame. Besides, the amplitudes squared are 
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In the case of the a s channel, the momentum transfer in the center of mass frame between 

the dark matter particles, g^, is given by 

1 V I 

^ — 2("^a "I" "^s ~ — m^i)^) + - cos^Y (4m^ + — m^)^ — 16m^m^ , (16) 

while in the a a and s s channels the corresponding q^ is given by setting m^ = ma in the 
above formula. Then, we can use the above formulas to get the velocity-averaged annihi- 
lation cross sections at freezeout or at present. For non-relativistic dark matter particles, 
we can expand the velocity-dependent terms by (cf )j ~ + hiv"^ and get the thermal cross 
sections as (o'f)j ~ + 6biTf/m^ with Tf ~ m^/20 being the freezeout temperature. We 
note however that, when the dark matter and scalar masses are almost degenerate and thus 
lead to a narrow box, the velocity expansion breaks down at freezeout and a thermal average 
of the full annihilation cross section is required. At zero temperature, we only have to keep 
the velocity-independent terms. 

In the non-relativistic limit, the velocity-averaged annihilation cross sections of the pro- 
cesses XX ~^ otf) 55 read 
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Note that, unlike XX ~^ '^■5, the annihilation channels XX ~^ (^(^^ present a strong phase 
space suppression when ma,ms — m^. This will be explicitly illustrated in Fig. [2] below. 

We include for completeness the s-channel dark matter annihilation cross sections into 
electroweak gauge bosons. The corresponding cross sections with axion mediator for a 
photon pair and Zj are given by 



1 2 2 16m^ 
{av),, = \c,,\ (4^2 _^2)2 + r2m2 ' (20) 

1 2| |2 IQm^ / m| \3 

^^^^^-^ - I'^-'l iAml-mir + TlmlV~^J ' ^^^^ 

which can be sizable when ~ 2m^. In the above expressions is the axion total decay 
width. The cross sections for the other s-channels with axion mediator ai^ 

_ 1 21 12 ^6m^ / m|\3/2 

^^''^^^ - 2^'^^"'^^' {Ami - <y + ^l^l ' ^ " ' ^^^^ 

1 ,2,„ ,2 IK A 



- T-|Ax| Ich/w^I —— — — ^T^— ^ 1-— ^ • (23) 



A-K {Am?, — m^y + m'^T'^ \ m 



There are in addition s-channel annihilation channels mediated by CP-even scalars. How- 
ever, as argued in Refs. [MIET], these are p-wave suppressed and thus irrelevant for indirect 
dark matter detection. These channels could though be important in determining the relic 
density at freeze-out when the CP-even scalar singlet has a sizable mixing with the Standard 
Model Higgs [36l[37]. 

In this paper we will be interested in the regime where < 2m^ and will assume a 
small mixing between the CP-even scalar and the Standard Model Higgs. Therefore, the 
relic density is essentially determined by the annihilation cross sections of the processes 



XX ~^ as,aa,ss, given by Eqs. (10)-(12). As will be shown in the next Section, depending 
on the dark matter coupling to the axion A;^ as well as on m^, ma and m^, it is possible to 
reproduce the measured relic density for a wide range of masses for reasonable values of the 
coupling A;^ and without any fine-tuning of parameters. 



4 Box-shaped gamma-ray signatures 

The annihilation channels xX ~^ and xX o,^, produce prominent gamma-ray features 
from the decay in flight of the axion into two photons or into one photon and one Z boson. 
Let us discuss first the energy spectrum produced in the annihilation xX ~^ (^s. The decay 
of the axion a — )■ 77 produces two photons with identical energy in the rest frame of the 
axion, E^^^-p = ma/2. However, in the galactic frame, where the dark matter particles move 
non-relativistically, the photon energy reads 

E!;'^^ = ^E^;^^p{l-VaCoser\ (24) 
^Thc typos due to a factor 8 in the annihilation cross sections in Refs. [36l[37j were corrected in Ref. [38] . 
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where 7 = 1/y 1 — f^, Va is the axion velocity in c units and 6 is the angle between the 
direction of the axion and the direction of the photon. Concretely, for the XX ~^ ^.s channel 
the velocity of the axion is 



Pa I m^ — m2\"2 



Ea M ^ ml + 4m2 ' ' ' ^^^^ 



which gives a photon energy in the galactic frame 



2 22 
Tn~, ( ^ m~ — m 




^^^ = ^(1 + ^^) I 1 -^(1 + ^^^) I . (26) 



Since the axion decays isotropically, the resulting energy spectrum presents a box-shape 
structure with the photon energy ranging from E^"'^ to E^"'\ where E^"'^ = ^Am^{l ± 

1 ~ ^2^2 ) and A = H-(m^ — m^)/(4m^). Analogously, the energy of the photon produced 
in the annihilation xX ~^ ^.s followed by the decay a — ?■ reads 



^{Zy) ^ 1^ f ) ( 1 + " , ' ) 1 - cos 



(27) 

producing also a box-shaped spectrum with photon energies ranging from Ei^^^ to E^^^\ 
where E*^^"* = ^Am^(^l — ^1 ± ^1 — -j^^^ ■ Therefore, the energy spectrum of hard 
photons produced in the annihilation channel xX ~^ o,s is 



dE^ ^(-^t) _ ^(7^) 



Q{E^ - E'y>)e{E\!^> - E^)BR{a 77) 



+ -T^:^y—T^QiE, - E^'^^mE^^'^ - E,)BR{a ^ Z7) , (28) 
E_^ — E_ 

where is the Heaviside function. In the presence of mixing between the singlet scalar and 
the Higgs boson, there are two annihilation channels: xX ~^ ^.s and xX ~^ (^h. In that case, 
the spectrum can be derived from the previous expressions substituting rris by the mass 
eigenvalues associated to s and h. For a small Higgs mixing, the dominant annihilation 
channel is xX ~^ '^5. In our analysis we assume that this is the case. 

When the dark matter particle annihilates in the channel xX ~^ (^0,1 the energy of the 
photon produced in the decay in flight of the axion can be straightforwardly calculated from 
Eqs. upon replacing by ma, namely 



rni I „ / ml 



for the case of the decays a — )■ 77 and 
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Figure 2: Contours of the branching ratios of the processes XX ~^ (^^i ^ function of 

the mass ratios ma/m^ and ms/m-^ (top), and as a function of 5a = {m^ — ma)/m^ and 5s = 
{my.—ms)/m^ (bottom) in order to highlight the region of the parameter space where ms,ma ~ m^. 
In the top plots the shaded grey region is kinematically inaccessible, while in the bottom plots 
the hatched region marks the parameter space where the width of the box (either from a — )• 77 or 
a Z7) is 10% or less of the maximum energy of the photons, i.e. — < Q.lEj^. We also 
show in both panels the location of our exemplary points with ma/rrij^ = mslra^ = 0.999 for the 
narrow box scenario ("A") and 0.1 for the wide box scenario (a). 



for the case of the decays a — )■ Z7. The energy spectrum of photons produced in the 



annihilations is then given by Eq. ( 28 ) with an additional factor of two due to the fact that 
two axions are produced in each annihilation. 

Finally, the annihilations of Dirac dark matter particles produce a gamma-ray flux [H |2] 

^=i^5:<->/^As^A/"^... . (31) 

where (erf) / is the velocity-averaged annihilation cross section in the channel /, Af2 is the 
observed field of view and Jann = Jios p'^m the integral of the squared dark matter 



density pdm along the line of sight. The convolution of the flux in Eq. (31) with the 
experimental energy resolution is done following Ref. [lO] (see details below). 

For a given annihilation channel, the energy spectrum of photons depends on ma/m^, 
ms/m^, BR(a — )■ 77) and BR(a — )■ Zj). The cross sections for the annihilations xX ~^ 
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as, aa, ss relative to the total cross section are shown in Fig. |2]for various values of ma/my^ 
and ms/m^. It is also shown the region of the parameter space where the width of the box 
(either from a — )■ 77 or a — )■ Z7) is 10% or less of the maximum energy of the photons, 
— E_ < O.IE^, which corresponds to the typical energy resolution of a gamma-ray 
telescope at the endpoint of the box. One can then identify two possible scenarios depending 
on the width of the box compared to the experimental energy resolution: the "narrow box" 
scenario and the "wide box" scenario. In the former case, the axion is produced almost 
at rest and the gamma-ray spectrum resembles a line, while in the latter the axion is very 
relativistic and the photons produced in its decay present a large spread in energy. In both 
cases the measured relic density can be reproduced for a wide range of dark matter masses 
with A;^ ~ C'(l), as demonstrated in Fig. |3} In Fig. |4| we also depict the ratio of the 
present to thermal total annihilation cross sections, {o'v)o/ {av)th, in the parameter space Ss 
vs 5a- We find that the narrower the gamma-ray box becomes, the more sizable the velocity- 
dependent contribution to the dark matter annihilation cross section and the smaller the 
ratio of the present to thermal cross sections are. For instance, we get {<7v)o/ {av)th — 0.13 
for a narrow box with ma/m^ = nis/m^ = 0.999 while {(Jv)o/ {<7v)th — 0.96 for a wide box 
with nia/m^ = nis/m^ = 0.1. 
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Figure 3: Contours of the total annihilation cross section at freeze-out in units of the thermal 
cross section for Dirac dark matter particles, (erf )th = 6 x 10~^^ cm'^ s~^, as a function of and 
A^. The left (right) panel shows the narrow (wide) box scenario with majm^ = ms/rriy. = 0.999 
(0.1). 

We have calculated limits on the annihilation cross section as a function of the dark 
matter mass for the process XX ~^ '^'^ in the narrow and wide box scenarios from the 
Fermi-LAT observations of the galactic centre region (Reg3 in and from the H.E.S.S. 
observations of the galactic ridge region [H] . The energy resolution of Fermi-LAT is taken 
from Refs. |l5lll6], while for H.E.S.S. a constant relative energy resolution of 15% is assumed. 
Besides, for the dark matter distribution we adopt the Einasto profile, favoured by recent 
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Figure 4: Contours of the ratio of the present to thermal total annihilation cross sections, as a 
function of 5a = {m^ — ma) /rrix and 5s = {m^ — ms^jm^. The hatched region marks the parameter 
space where the width of the box is 10% or less than the maximum energy of the photons. We 
also show the location of our exemplary points with ma/m^ = ms/m^ = 0.999 for the narrow box 
scenario (-^) and 0.1 for the wide box scenario (a). 

A^-body simulations [SI SHI SS] , 

(32) 

with Ts = 20 kpc and a = 0.17 [501 El] and normalised to pBui^ = 8.5 kpc) = 0.4 GeV cm~^ 
[S21 ESI El ES]. The J-factors read ^ /^^ dQ Jann ^ 3.02 x 10^^ GeV^ cm'^ for the galactic 
centre region (Reg3 in [TT]) and ^ j^^dVlJs,rm — 1-29 x 10^^ GeV^ cm"^ for the galactic 
ridge region defined in [H] . The 2a limits are shown in Fig. [sj assuming for concreteness 
ma/frt^ = nis/m^ = 0.999 in the case of the narrow box scenario and ma/m^ = m^/m^ = 0.1 
for the wide box scenario (choosing a different mass for a and s would lead to very similar 
phenomenology). These values for the mass ratios are indicated in Figs. |2] and |4] with a 
star and a triangle, respectively. Furthermore, in Fig. |5] we assume BR(a — )■ 77) = 1, 
which holds to a good approximation when < mw, as commonly occurs in the wide box 
scenario, but the limits can be trivially rescaled for a different branching ratio. We have 
computed the limits following two different approaches to bracket the uncertainty regarding 
background modelling: the solid curves labelled "conservative" correspond to the limits 
assuming that the background flux is zero, while the dashed curves labelled "aggressive" 
assume that the background meets the central points of the measurement. Let us note that 
these limits apply not only to the model discussed in this paper but to any model where the 
dark matter particle is a Dirac fermion that annihilates into two (pseudo-)scalars, which in 
turn decay in flight into two photons. In the case of a Majorana dark matter particle, the 
limits in Fig. Is] are stronger by a factor of two (and the thermal cross section is reduced by 
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Figure 5: Upper limits on the annihilation cross section in the channel XX ~^ (^(^ from Fermi- 
LAT observations of the galactic centre region |17| and from H.E.S.S. observations of the galactic 
ridge region [H], assuming BR(a — )• 77) = 1. The left (right) panel shows the narrow (wide) box 
scenario with ma/m^ = ms/m^ = 0.999 (0.1). The solid and dashed curves indicate, respectively, 
the limits following the conservative and aggressive approaches described in the text. The thermal 
cross section for Dirac dark matter particles, (cjf )th = 6 x 10~^^ cm^ s~^, is shown by the horizontal 
line. 

a factor of two). The limits corresponding to the decay of the (pseudo-)scalars into Z7 can 
be readily derived from these plots with the appropriate substitutions in the kinematics, as 
well as the limits for the annihilation channel xX ~^ ^.s. 

We illustrate the potential of axion-mediated dark matter to produce observable sig- 
natures by showing the predicted gamma-ray flux for two exemplary points, one for the 
narrow box scenario and one for the wide box scenario, which produce an excess at around 
130 GeV, motivated by the recent hint of a sharp feature at this energy [121 [H] . Although it 
is tempting to try to explain this signal in the framework of our model, we will not attempt 
to make a fit to the data. 

We consider first a narrow box scenario corresponding to = 250 GeV, = = 
0.999mjj, = 249.75 GeV, shown as a star in Figs. [2]and|4| and which gives BR(xx — )■ as) ~ 
1.0, BR(xx — )■ aa) ~ 0, and {(Jv)o/ {av)th — 0.13. Furthermore, we consider the case with 
Ci/c2 = 3, such that BR(a — )■ 77) ~ 0.4. The total flux predicted by the model for the 
region of the sky Reg3 defined in [T7j with this choice of parameters is shown in Fig. |6] (left) 
as a thick red line. Here, we adopted an Einasto profile and a background flux following 
the simple power law = 1.5 x 10~^(i?^/GeV)"^-^^ GeV~^ cm"^ s"^ sr~^, which fits well 
the data below 70 GeV. The model predicts a gamma-ray feature that stands out over the 
power-law background at ~130 GeV without invoking any boost factor. For other values of 
Ci/c2 the gamma- ray feature can become fainter or stronger, however this change in intensity 
could be compensated by moderate variations of the astrophysical parameters, e.g. a different 
choice of halo profile or an astrophysical boost factor that could enhance the signal by a 
factor of 0{1). Note that for annihilations xX as a spectrum of continuum photons can 
be generated by the decays of s and a, and therefore the corresponding constraints should 
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Figure 6: Predicted gamma-ray spectrum for a narrow box scenario with rriy. = 250 GeV, = 
uia = 249.75 GeV assuming ci/c2 = 3 (left) and for a wide box scenario with = 150 GeV, 
TTia = rus = 15 GeV (right). In both panels, the dark matter coupling A,^ is determined as a 
function of dark matter mass by the thermal annihilation cross section {av)th = 6 x 10~^^ cm^ 
(see Fig. [s]) and the present annihilation cross section follows from Fig. |4j The branching ratio for 
the different annihilation channels and the branching ratios for the axion decay are predicted by 
the model and can be read from Figs. [2] and [T| respectively. The background flux is assumed to be 
a power law in the energy range shown and has been chosen to fit the data below 70 GeV (black 
line) . 



be taken into account. These constraints on continuum photons typically place an upper 
limit of (9(10 — 100) in the ratio of total cross section to line cross section |56l |571 158], so 
the models considered here are safe as long as the branching ratio of the process a — )■ 77 is 
sizable. 

The narrow box case analysed in the previous paragraph has fairly degenerate dark 
matter and axion masses and could arguably be considered fine-tuned. Nevertheless, the 
expectation of a gamma-ray excess holds in a much larger region of the parameter space. 
We illustrate this fact using a second exemplary choice of parameters, corresponding to 
a wide box scenario and which in fact produces a gamma-ray flux overshooting the data. 
Concretely, we have taken = 150 GeV, nia = rus = O.lm^ = 15 GeV, such that the 
axion decays dominantly into two photons, i.e. BR(a — )■ 77) ~ 1. This point is shown as 
a triangle in Figs. [2] and |4] and yields BR(xx — )■ as) ~ 0.25, BR(xx ~^ Q;a) ~ 0.25 and 
{crv)o/ {av)th — 0.96. The total gamma-ray flux that results from this choice of parameters 
is shown in Fig. [6] (right) as a thick red line and turns out to be too intense to reproduce 
the Fermi-LAT observations. Again, moderate variations of the astrophysical parameters 
could lower the intensity of the signal producing a gamma-ray flux in better agreement with 
the data. However, and as we already mentioned above, we will not attempt to explain 
the 130 GeV excess within this framework. We merely present this example to illustrate 
that models producing gamma-ray boxes can generate, for the thermal cross section and 
typical astrophysical parameters, spectral features intense enough to be probed in present 
and future gamma-ray telescopes. 
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5 Conclusions 



With several instruments providing observations of exquisite quality, it is now possible to 
search for fine spectral features in the gamma-ray emission from our Galaxy and beyond. 
This opportunity - not present several years back - is of extreme importance for dark matter 
searches and is perhaps one of the cleanest ways to glean on dark matter phenomenology. In 
the present contribution, we have focused on a particular type of spectral features, gamma- 
ray boxes, which stand out over the background at high energies. Interestingly, gamma-ray 
boxes are naturally realised in dark matter models without the need for fine-tuned parame- 
ters. Constructing a complete framework for axion-mediated dark matter, we show explicitly 
that both narrow and wide boxes are possible configurations of the photon spectrum induced 
by annihilations. Either way, thermal relics unambiguously foresee observable gamma-ray 
fluxes. On a more general note, we have also obtained model-independent constraints on 
box-shaped features from dark matter annihilation using the data provided by Fermi-LAT 
and H.E.S.S.. Our results - that are, to the best of our knowledge, the first of their kind 
to extend to the multi-TeV range - are remarkably strong despite our straightforward data 
analysis procedure. A dedicated box search would be an important step towards exploring 
the full potential of the current gamma-ray observations and would probably improve upon 
the constraints presented here. More importantly, judging from our findings, the prospects 
for future lACTs look very exciting given that such instruments will likely probe thermal 
cross sections up to several tens of TeV. The search for gamma-ray boxes in lACTs (in 
particular, CTA) rests thus as one of the last hopes to test thermal relics at high masses, a 
regime where other dark matter signals become increasingly dim, rendering direct detection, 
collider searches and other indirect detection channels of little relevance. 
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